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ABSTRACT 

Aims. The objective of this study was to examine the kinematics of coronal mass ejections (CMEs) using EUV and coronagraph 
images, and to make a quantitative comparison with a number of theoretical models. One particular aim was to investigate the accel- 
eration profile of CMEs in the low corona. 

Methods. We selected two CME events for this study, which occurred on 2006 December 17 (CME06) and 2007 December 31 
(CME07). CME06 was observed using the EIT and LASCO instruments on-board SOHO, while CME07 was observed using the 
SECCHI imaging suite on STEREO. The first step of the analysis was to track the motion of each CME front and derive its velocity 
and acceleration. We then compared the observational kinematics, along with the information of the associated X-ray emissions from 
GOES and RHESSI, with the kinematics proposed by three CME models (catastrophe, breakout and toroidal instability). 
Results. We found that CME06 lasted over eight hours while CME07 released its energy in less than three hours. After the erup- 
tion, both CMEs were briefly slowed down before being accelerated again. The peak accelerations during the re- acceleration phase 
coincided with the peak soft X-ray emissions for both CMEs. Their values were ~60 m s~^ for CME06 and ~600 m s~^ for CME07. 
CME07 reached a maximum speed of over 1000 km s~^ before being slowed down to propagate away at a constant, final speed of 
~700 km s~^ CME06 did not reach a constant speed but was moving at a small acceleration by the end of the observation. Our 
comparison with the theories suggested that CME06 can be best described by a hybrid of the catastrophe model and breakout model 
while the characteristics of CME07 were most consistent with the breakout model. Based on the catastrophe model, we deduced that 
the reconnection rate in the current sheet for CME06 was intermediate, the onset of its eruption occurred at a height of ~200 Mm, and 
the Alfven speed and the magnetic field strength at this height were approximately 130-250 km s~^ and 7 Gauss, respectively. 

Key words. Sun: coronal mass ejections (CMEs) - Sun: flares - Sun: magnetic fields - Sun: corona - Sun: atmosphere 
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1. Introduction 

Coronal mass ejections (CMEs) are the ejections of large amount 
of mass and magnetic flux from the Sun to interplanetary space. 
The energy released during the process is of the order of 
, 10^^ - 10^^ erg. Statistical studies have reported variations in 
■ the observed propertie s and projected kinem atics of CMEs (e.g., 
' ICremades & Bothmeii 120041: IZhang et al.L ^004) Nevertheless, 
CMEs often rise with a small initial speed of a few kilome- 
ters per second, followed by a rapid upward expansion reaching 
several hundred or several thousand kilometers per second, and, 
eventually, propagate through the interplanetary space at the am- 
bient solar- wind speed. 

Several models have been proposed to explain the driving 
mechanism and observed properties of the CMEs. A valid model 
must be able to produce the observed kinematics, dynamics and 
properties of CMEs. In Sec. [21 we give a review of three repre- 
sentative models: breakout model (BO) ( Antiochos et al., 1999), 
catastrophe model ( CA) ( see, e.g., van Tend & Kuperus, 1978t 
Forbes & Isenbergi 119911: Forbes & Priest .1995: Lin et al 



and toroidal ins tability model (TI) (see,'e.g.,|Chen, 1989; 
Kliem & Torok, 2006), which have been investigated in this 
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study. Although the profile of the CME kinematics predicted by 
each model is qualitatively consistent with the observed three- 
phase profile (i.e., slow rise, eruption and slow down), the ex- 
act profile of each model is diflTerent. Several simulations also 
showed that a same model can produce diflTerent kinematic pro- 
files by s imply varying the initial conditions and magnet i c env i- 
ronment (iLvnch et al.Ll2QQ8l:ISchriiver et al.Ll2QQ8l:ICheil 119891) . 

In addition to the discrepancies among theoretical pre- 
dictions, observational studies that were not based on any 



model, also reported diflferent functi onal forms tha t best-fi tted 
their respective data. For instance, 'Sheele y et al.l (Il999h fit- 
ted their data by r(t) = ro + 2ro In cosh [v a(t -\- Ato)/2r a]', 
iGallagher et al.l (120031) formulated a double exponential func- 
tion a(t) = ^a~^ exp(-t/Tr) + a'^^ exp(^/Tj)j to describe the fast 
r ising and decaying acceleration profile; and [Alexander et al.l 
(I2002h demonstrated that their fast CME was best fitted by the 
polynomial, h(t) = ho -\- vot + ct^, with m ^ 3 .7, which is con- 
sistent with the results by lSchrijver et al.l (l2008h . Although these 
studies were able to derive how the velocity and acceleration 
evolve through the duration of CMEs, they did not provide a 
physical explanation and driving mechanism behind such kine- 
matic profiles. Besides, while most of these studies chose poly- 
nomials or exponential functions or a simple combination of the 
two, there is no guarantee that the chosen function is the only 
one that can fit the observational profile. 

Studies to verify a specific model by comparing observations 
with models have also been carried out. For example, by select- 
ing o nly th o se tha t can be characterized as flux-rope type CMEs, 
Krall et af] (l200ll) showed that the kinematic curves predicted 
by the TI model matched the kinematics in the higher corona 
but not the low er coronal region. An extensive examination of 
a fast CME by iManoharan & Kundul (l2003h based on the data 
from various instruments showed supportive evidence for the 
BO model. However, Bong et al. (2006) reported that while the 
magnetic-field configuration and X-ray loop field connectivity in 
their data indicated a breakout process, the acceleration profile 
they obtained contradic ted the BO model. The studies to ver- 
ify the CA model (e.g., 'Alexan der et al.L l2002t ISchrijver et al.L 
.2008.) . however, often deduced their conclusions by compar- 
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ing the early-stage expression of CA model with the accelera- 
tion phase and/or the late stage of the observed CME kinematic 
curves. As already been demonstrated in all observational stud- 
ies, the kinematics of CMEs can change greatly through the du- 
ration of the event. Hence, it is unsurprising that the functional 
form for the early stage does not match with the profiles of the 
later stages. 

The objectives of this paper are first to examine the features 
and kinematics of CME events, and then to determine the pos- 
sible mechanism behind the observed phenomena by comparing 
the observations with various theories. The theory most consis- 
tent with the observation can thus provide information on the 
magnetic environment where the eruption occurs. We have cho- 
sen two CME events for this study. We first derived their kine- 
matic profiles and examined the associated X-ray emissions and 
flares, and then compared the results with those derived from 
the three difl'erent models. The kinematic profiles from the mod- 
els have been mainly obtained from simulations. The only stage 
of which an analytical expression of the kinematics is available 
is the early stage of the eruption. Hence, we first focused on a 
qualitative comparison of the characteristics between the obser- 
vations and diff'erent models, and then carried out an additional 
quantitative examination by fitting the observational data with 
a number of model expressions. The objectives of the quantita- 
tive examination were to verify the assumptions employed in the 
models to derive these expressions, and/or to obtain an empirical 
expression for the observed kinematics. The values of the fitting 
parameters may be used to infer the magnetic environment and 
the initial conditions at the onset of the eruption. 

The rest of the paper is organized as follows: the three mod- 
els are reviewed in Sec.O The two selected CME events and how 
their motions were traced are described in Sec.O The qualitative 
and quantitative comparisons with the models are explained and 
discussed in Secs. l3.3l and l3.4[ We conclude our study in Sec.O 



2. Review of CME models 

2.1. Catastrophe model (CA) 

In a two-dim ensional CA model (e.g, iForbes & Isenbergi Il99ll: 
iPriest & Forb es, 2002), the flux rope is represented by a ring cir- 
cling the Sun, and, hence, does not have anchoring ends. This 
flux-rope magnetic-field configuration is generated by placing 
two line sources with opposite polarities on the surface around 
the Sun. The magnetic field line connecting the two sources is 
initially balanced by a downward force from the tension of the 
field line and a upward force due to the magnetic pressure gradi- 
ent. If the photospheric motion moves the two sources closer, 
the tension force would increase, which pulls the field lines 
downward. However, the downward movement, compressing the 
field lines, subsequently increases the pressure force. When the 
sources move closer than a critical distance, the pressure force 
takes over and the flux rope is propelled and accelerated upward, 
resulting in the eruption. The model shows that, at the critical 
point, the distance of the sources {d) and the height of the flux 
rope Qi) is related hy h = d/2 = Aq. After the rope rises to a cer- 
tain height (h ^ d ^ 2Ao), a current sheet begins to form below 
it, which acts as a drag force against the rising of the rope. Using 
their 2-D model implemented with an isothermal atmosphere, 
iLin & Forbesl (l2000l) showed that the kinematics of CMEs is de- 
pendent on the reconnection rate in the current sheet. In their 
study, the reconnection rate is constant along the current sheet, 
and is prescribed by the inflow Alfven Mach number. Ma, at the 
midpoint of the current sheet. They found two critical values of 



Ma, 0.005 and 0.041, dividing CMEs into three classes: oscil- 
latory, with and without a deceleration phase. Since oscillatory 
CMEs have never been observed, the first critical value essen- 
tially provides a lower limit of Ma for a CME to happen. If the 
reconnection rate is between the two critical values, the CME 
in their simulation first undergoes a fast acceleration phase for 
~12 min after the loss of mechanical equilibrium, followed by a 
deceleration period = 20 - 100 min) and then is re-accelerated 
again (cf. Fig. 6 in the paper). For CMEs with a high reconnec- 
tion rate (Ma > 0.041), the deceleration period would not occur. 
To examine the temporal behaviour of the reco nnection rate and 
the rate of magnetic energy release ( dWry^c, a /dt), \FnQst & Forbesl 
( 2000) used the electric field at the X point, Eq, to represent the 
reconnection rate, and allowed it to freely evolve. We can see 
from their result (Fig. 11.6 in Priest & F orbes (2000)) that the 
profile of Eq closely follows the profile of dM^^^jdt and that 
the peaks of the two profiles almost coincide. In other words, 
we may infer from their result that the peak of the reconnection 
rate occur at a similar time as the peak of dW^^^gldt. Based on 
the above results and assuming that the rate of thermal-energy 
change {dWi\^ldt) follow s the rate of magnetic-energy release 
(iReeves & Forbesl [2005 ") and that dWi\^ldt can be represented by 
the soft X-ray (SXR) emission light curve, we can deduce that 
the SXR may reach its emission peak when the reconnection rate 
is the highest. We can also deduce from the CA theory that the 
SXR and hard X-ray (HXR) emissions, which are indications 
of reconnections and flares, should begin to increase/occur later 
than the launching of the CME (i.e., after the reconnections start 
in the current sheet). 

An analytical expression for the kinematics of a thin flux 
rope (i.e., the radius of the flux rope a ^ Qi) before the forma- 
tion of the cu rrent she et (i.e, hlA{) ^ 2) c an be derived (e.g., 
iLin & Forbesl 2000: Priest & Forbesl l2000h : 



h 



ln(A] + ^-2tan-^fA 



1/2 



(1) 



where h is the height of the flux rope, and ho and ho correspond to 
the initial height and speed, respectively. Ao is the critical height 
(i.e., the height at the critical point), and Vaq is th e Alfven speed 
at /z = Ao. Starting with Eq.[TJ Prie st & Forbesl (^2000) further 
derived simplified expressions for the "early" and "late" part of 
this stage. By considering the "early" part as before the time 
reaches the Alfven time scale, i.e., t <^ Ao/Vao, they obtain: 



h ^ Ao -\- hot + 



4y. 
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The expression for the "late" part was derived by assuming 
h/Ao » 1 but I lnh\ still less than | lna\\ 



Vao 



In 



1/2 



(3) 



After the formation of the current sheet, the system becomes 
too complicated for an analytical expression of the CME kine- 
matics to be derived. 



2.2. Toroidal instability model (Tl) 

The TI model dCtol fT989l: iKliem & TorSfl [2QQ6h considers a 
pre-existing coronal flux rope that is initially in a stable equilib- 
rium, and aims to explain what initiates the eruption of the loop. 
There are two dominant forces acting on the flux rope. One is 
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an outward Lorentz self-force (hoop force) resulting from the 
current flowing in the flux rope, and the other is a compressing 
Lorentz force due to the interaction between the poloidal mag- 
netic fields outside the flux rope and the current inside the flux 
rope. If the external magnetic field (^ex) decreases sufficiently 
fast from the flux rope and vanishes in the infinity, the compress- 
ing Lorentz force would not be strong enough to counter-act a 
sudden increase of the hoop force resulting from, for instance, a 
surge of the current in the flux rope. The flux rope would then 
rapidly and continuously expand, leading to a CME eruption. 
Th e condition of ^ ex for the instability to happen was derived 
bv lBatemanI (119781) : -RdlnB^^/dR > 3/2, where 7? is the major 
radius of the flux rope. We can infer from the above descrip- 
tion that TI theory does not require reconnections to accelerate 
a CME, which is a major distinction between TI and the other 
two models. Instead, TI proposes that a flux rope can erupt and 
become CME as long as the external magnetic fields satisfy the 
instability condition. 

The equation of motion derived from the TI model is too 
comp licated to be solved analytically (see, e.g. jKliem & Torofl 
l2006h . However, by considering only two counter- acting Lorentz 
forces on the rope, ignoring temporal changes in the external 
field and its resulting flux in the rope, and assu ming a simple 
profile for the external field (B^^H) = M""), Kliem & Torokl 
E0Q6) showed that the expression for the beginning stage of the 
instability can be approximated as a hyperbolic function: 
p 

h(T) = y sinh(PiT), h = H/Ho - 1 « 1, (4) 

where H and Hq are the height of the rope and the rope height at 
the on- set of instability; r is the time normalized by Alfven time; 
Po is a parameter composed of velocity, inductance, rope height 
and rope radius at the on- set of instability; and Pi = ^Jn - ricr 
is associated with the profile of the external magnetic field. 
iKliem & Torokl (l2QQ6 l) also presented the profiles from their nu- 
merical simulations based on the three aforementioned simpli- 
fications (i.e., only two acting forces, a constant external field, 
and a simple profile for the external field). One distinctive fea- 
ture in their simulation results is that the acceler ation all shows 
a fast rise and a more gradual decay. However. ISchrijver et al.l 
(|2QQ8|) demonstrated that the height-time profile of a TI model 
can change from a hyperbolic function to polynomials by simply 
tuning the initial conditions, which would change the accelera- 
tion profile from a fast initial rise to a more gradual one. 

2.3. Breakout model (BO) 

One d istinctive requirement for the BO model (lAntiochos et al.L 
fT999h is that the flux rope is formed in a multi-polar arcade sys- 
tem, which consists of a low-lying central arcade, two low-lying 
side arcades (one on each side of the central arcade), and a large 
arcade overlying the three arcades. The central arcade is what 
would later become the CME flux rope. The shearing of the cen- 
tral arcade boundaries causes the arcade to rise, which subse- 
quently results in the reconnections with the overlying arcade. 
As this reconnection (the breakout reconnection) removes some 
of the overlying magnetic field lines, the rising of central arcade 
becomes faster, and, in turn, triggers more and faster breakout 
reconnections. This process manifests as the eruption phase in a 
BO scenario. As the flux rope rises, it draws the magnetic field 
lines of the central arcade close to each other to form a current 
sheet underneath. When the current sheet reaches a certain crit- 
ical dimension, tearing instability would cause a second set of 
reconnections, which not only release and further accelerate the 



CME rope, but also result in post-CME flares. After this set of 
reconnections cut ofl' the current sheet, the magnetic field lines 
of the side arcades would move in to form another current sheet, 
in which a third set of reconnections would happen to reform 
and restore the magnetic fields after the eruption and also trigger 
another series of flaring activities, which they identified as the 
ribbon flares. 

Since the breakout reconnection is what allows the flux rope 
to rise and erupt in a BO picture, we expect that, ideally, signifi- 
cant SXR and HXR emissions should be detected from the very 
beginning of the loop rising. If the second set of reconnections, 
responsible for cutting off' the tethering force from the current 
sheet below the CME, can be equated to the reconnection pro- 
cess in the CA model, the temporal behaviour of the reconnec- 
tion rate during this phase may follow the proffle of the corre- 
sponding SXR emissions, as we discussed in Sec. 12.11 Because 
of the reconnections with both the overlying ffeld lines and the 
trailing current sheets, we can expect indications of ffaring activ- 
ities to be detected both below and on the sides of the CME loop 
during diff'erent stages of the event. However, we must keep in 
mind that the emissions and ffares during some stages might not 
be detectable due to observational limitations. 

Although the BO model does not have an analytical ex- 
pression for the kinematics either, there are t hree possib l e pro- 
files produced by three diff'er ent simulations dLvnch et a l.ll2004l 
I2QQ8I: iDeVore & AntiochosL l2QQ8l) . While iLvnch et a l. (2004) 
showed that their 2.5 dimensional simulation result can be fit- 
ted by a single, constant- acceleration profile: 

1 9 

h(t) = /zo + vo(t - to) + -ao(t - to)\ (5) 

the three-dimensional simulation by iLynch et aP (l2008l) pro- 
duced a more complex kinematic profile with three phases of 
constant acceleration: a long interval of low, constant accelera- 
tion during the foot-point shearing and breakout reconnection, 
followed by a short interval of high, constant acceleration dur- 
ing the second set of reconnections (i.e., reconnections in the 
current sheet formed by the central-arcade magnetic field lines), 
and finished with a long inter val of constant velocity which sim- 
ply means zero acceleration. iDeVore & Antiochoi (l2008l) pro- 
duced a slightly diff'erent three-phase acceleration profile from 
their simulations: a long interval of slow acceleration during the 
shearing of the footpoints, a short interval of fast acceleration 
through the breakout and ffare reconnections, and a short inter- 
val of fast deceleration at the end indicating the restoration and 
reformation of the magnetic ffeld conffguration. 

3. Observations and data analysis 

The two CMEs selected for this study were observed on 2006 
December 17 (CME06) and 2007 December 31 (CME07). Both 
events were initiated at the solar limb, which allow their radial 
motion to be easily observed. 

We have utilized observations f rom EIT (E xtreme ultravi- 
olet Imaging Telescope; Delaboudiniere et al.l ([l995)) 195 A 
wavelength ba nd and LASCO (L a rge A ngle and Spectrometric 
Coronagraph; iBrueckner et al.l (Il995h ) for CME06, and 
STEREO B (Solar TErrestrial RElation Observatory; 
How ard et all (l2008l) ) instruments (EUVI 171 A wavelength 
band, corl and cor2) for CME07. The X-ray emission data were 
obtained from GO ES (G eostationary Operational Environmental 
Satellite; iGarcial (11994 )) and RHESSI (Reuven Ramaty High 
Energy Solar Spectroscopic Imager: . Lin et al.. (200Z) ). 
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Fig. 1. The image sequence of CME06 observed by EIT and 
LASCO C2 and C3. The EIT images have been rotated so that 
the loop top moves along X axis, as indicated by the dash line. 
The pivot of the rotation is plotted as a star, which is located at 
X~140 along the dash line. The crosses in all images mark our 
determined CME front. 



We first applied a temporal running diff'erence to the time- 
series images to enhance the contrast, and then determined the 
height vs. time relation of the CME by tracing the motion of 
the CME front. Our determined CME fronts of the two events 
are plotted as crosses in Fig. [T] for CME06 and Fig. [3] and |4] for 
CME07. Since the accuracy of the height determination is the 
most important factor for the accuracy of the rest of the analysis, 
we repeated the height determination at least ten times, and cal- 
culated the standard deviation of the multiple trials. The standard 
deviations were of the order of lO^km for EIT and EUVI, and 
10^ - lO^km for LASCO and STEREO corl and cor2. We then 
used the standard deviation plus the spatial resolution of the im- 
ages as the uncertainties. The observational velocity and accel- 
eration were subsequently calculated by taking the time deriva- 
tives of the height. The errors in the derived velocity and accel- 
eration were propagated from the height uncertainties using the 
IDL function derivsigC). The kinematic profiles of the two 
CMEs, along with the error bars and the GOES SXR profiles, 
are plotted in Fig. |2] (CME06) and Fig. [5] (CME07). In the left 
panels of both figures, we plotted only the results of the lower 
corona to allow better visibility of the kinematics in this region. 



3.1. 2006 December 17 CME (CME06) 

This CME was highly inclined in the low corona within the field 
of view of EIT. Hence, we rotated the EIT images with respect to 
the "base point". The base point is the center between the loop 
footpoints at the limb. The line from the loop top to the base 



point is called "base line". This rotation resulted in the apparent 
propagation of the loop top along the X-axis. Three examples of 
the rotated EIT running-diff'erence images are shown in the top 
row of Fig. [U The base line is plotted as a dashed line, and the 
base point is marked as a star. The images of LASCO are not 
rotated because the CME appears to move in the radial direction 
at such height. 

The CME, initially observed by EIT, was embedded among 
other loops until it began to rise shortly after 12:00 UT. 
The CME was within the field of view (FOV) of EIT until 
-15:12 UT. It was later observed by LASCO C2 between 15:30 
and 19:00 UT and by C3 between 17:00 and 21:18 UT, after 
which it became too faint and diff'use to be reliably identified. 

The derived kinematics of the CME (Fig. [2]) shows that the 
rise of the CME was initially at a constant speed of a few km s"^ 
until ~ 14:30 UT, when the CME suddenly became eruptive, indi- 
cated by the steep rise of velocity and acceleration. The eruption 
was followed by a quick deceleration, reaching ~ -10 m s"^. 
Although the exact duration of the eruptive phase is uncertain 
due to the gap between the EIT and LASCO C2 FOV, the ac- 
celeration profile in Fig. [2] indicates that the CME turned from 
eruption to deceleration in ~1.5 hr. After the deceleration, the 
CME was accelerated again, and then slowed down to propagate 
away at a small constant acceleration. 

The GOES SXR emissions exhibited two rising phases, an 
initial steep rise followed by a more gradual rise. The first rising 
phase began when the CME almost reached the boundary of the 
EIT FOV (> 14:50 UT). The second rising phase, as can be seen 
in Fig. [21 is where the first CME acceleration peak occurred. The 
SXR reached its peak at ~ 17: 15 UT, when the CME reached the 
re-acceleration peak and a flaring loop was detected by EIT. We 
also examined the relevant Hinode XRT images, which showed 
that this flaring loop was located directly below the erupted loop 
and that it contracted slightly before rising to form a cusp- shape 
loop top. By comparing the RHESSI and the EIT images, we 
found that the sources of the RHESSI emission (3-6 keV and 6- 
12 keV) were located at the top of this post-eruption flaring loop. 
The GOES SXR emission finally decreased to a background 
magnitude after 21:00 UT. 

3.2. 2007 December 31 CME (CME07) 

Fig.[3]shows the early stages of the CME emergence as observed 
by EUVI. The figure revealed that a system of loop arcades ap- 
peared at -00:41 UT, and quickly rose radially until 00:48 UT, 
after which the top of the arcades became visibly connected to 
form the front of a balloon- shaped CME, with a round top and a 
narrow lower part connected to the solar surface. EUVI also cap- 
tured some flaring activities and foot-point brightenings a few 
minutes before the emergence of the CME. After the emergence. 
Fig. |4] shows that the CME began to quickly expand sideways 
while the radial rise became slower than previously. The CME 
retained the round- shaped top throughout the EUVI and corl 
FOV, but changed to a shape resembling conjoined balloons by 
the time it entered the cor2 FOV. This suggests that the flux rope 
may have been turned, by solar rotation and/or footpoint shear- 
ing, to also show part of the farther side of the rope. 

The velocity profile of CME07 (second row in Fig. [5]) shows 
that the first reliable velocity point of the CME is ~200 km s"^ . 
Unlike CME06, the velocity of the CME07 did not show a clear 
drop, but rose monotonically to a peak value of over 1000 km s"^ 
in under 30 min. After reaching the peak, the velocity dropped to 
a final constant value of -700 km s"^ . Despite the monotonic rise 
of the velocity, the derived acceleration (bottom row of Fig. O 
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Fig. 2. The kinematic profiles of CME06. Diff'erent symbols rep- 
resent the observational results from diff'erent instruments, as de- 
noted in the upper right panel. The GOES X-ray emission pro- 
file has been scaled to enhance the viewing, and hence the exact 
magnitude is irrelevant. The right column shows the combined 
profile of all three instruments. The left column shows only the 
results from the low corona to give a better view of the kinemat- 
ics in this region. 



exhibited a decrease during the first 15 min, and then rose to a 
plateau value of ~600 m s"^, which lasted for ~15 min. After 

that, the acceleration dropped to 150 m s"^ before returning 

to almost zero at the end. We infer from the high initial velocity 
and the decrease of the acceleration at the beginning that the ob- 
servations did not detect the initial eruption phase, during which 
the CME may have reached an acceleration peak. 

The SXR emission from GOES began to rise shortly af- 
ter 00:40 UT, reached a peak around 01:10 UT, at which time 
the CME reached its peak velocity, and dropped to background 
magnitude roughly after 03:00 UT when the acceleration of the 
C ME reduced to almo st zero. The analysis of the same event 
bv lRafterv et"aD (l2010h . which utiHzed data from MESSENGER 
SAX and RHESSI, showed that the flare/SXR emissions actu- 
ally began much earlier at 00:27 UT and lasted until 03:50 UT. 
This indicates that the GOES observation was occulted during 
the beginning stage o f this event. A study o f the RHESSI HXR 
data of this event by iKrucker et al.l (l2009allbh revealed that the 
RHESSI HXR sources were located above the SXR sources and 
were detected at 00:48 UT and 01:02 UT. The authors deduced 
that the HXR sources may indicate the location of particle ac- 
celeration. Comparing the timing of these two HXR detections 
with the kinematics profile of this CME, we found that they oc- 
curred respectively at the time of the loop emergence and at the 
beginning of the re- acceleration. 



Fig. 3. The image sequence to show the formation stage of 
CME07 as observed by STEREO B EUVI 171 A. The crosses 
mark our determined CME front. 
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Fig. 4. The image sequence of CME07 observed by STEREO 
B instruments, EUVI 171 A, corl and cor2. The crosses in all 
images mark our determined CME front. 
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Fig. 5. The kinematic profiles of CME07. See the caption of 
Fig.[2]for the description of the figure. 



3.3. Qualitative examination 

The first step of this examination was to identify the properties 
that can be deduced/derived from and can be used to distinguish 
diff'erent models. These properties must also be derivable and/or 
detectable from the observation data we have. After these dis- 
tinctive properties were identified, we then examined the consis- 
tency between the features from our observations and the prop- 
erties of each models. 

A CME can be qualitatively identified as CA type if it shows 
the following features: the CME loop begins to rise before an 
increase in the X-ray emission is observed; CME- associated 
flares occur only after the launch of the CME; and the kine- 
matics r esemble any profile produced by the CA model (e.g., 
iLin & F orbes, 2000). If there is no flaring activities around the 
CME launch site and/or if the acceleration profile of a CME 
resembles any of the TI simulation results, it is likely to be a 
TI type. However, we would need a quantitative examination 
(i.e., fitting) to verify whether the profiles proposed by TI indeed 
match our observational results. To classify a CME as a BO type, 
we look for the following signs: significant X-ray emissions and 
flaring activities are detected from the very beginning of a CME 
event; there are magnetic arcades at the vicinity where the loop 
is formed; there are flare ribbons after the CME takes ofl'; and/or 
the kinematic profile resembles any of the simulation results. 

3.4. Quantitative Examination 

The mathematical profiles are often derived by imposing some 
approximations and assumptions. Changing the initial condi- 
tions and assumptions, or using diff'erent approximations, can 
often lead to a diff'erent theoretical profile without changing the 
model itself (e.g.. .Schrijver et al.».2008.) . Besides, the analytical 



expressions derivable from the theories are generally limited to 
the early stage of a CME process. An expression for the entire 
process of CME has not been available due to the complexity 
of the realistic models. Hence, the purpose of this quantitative 
examination is to verify the mathematical profiles derived from 
diff'erent CME models and, when the theoretical expression is 
unavailable, to empirically find a functional form for the kine- 
matics of the CME. Our strategy was to first apply several model 
profiles to fit the height-time data, then calculate the theoreti- 
cal velocity and acceleration by taking the derivatives of the fit- 
ting results, and lastly compare the theoretical and observational 
kinematic profiles. 

We implemented four fitting models for this exercise, CA, 
CAlate.e, TI and Poly fits. CA fit was of the form of Eq.[T] and 
the fitting parameters were the critical height Aq, initial speed 
vo, and Alfven speed at the critical height Vao- The profile for 
CAlate.e fit was Eq.[3]with a variable exponent, that is. 



h 



-Vao 

71 



ln|A)_^ 



(6) 



The fitting parameters for CAlate.e fit were the critical height 
/lo, Alfven speed at the critical height Vao and the exponent a. 
We also tested the original Eq.[3]by setting a equal to 0.5, and 
the result will be called CAlate. Because CA and CAlate.e fits 
were based on non-linear diff'erential equations, we conducted 
the fittings by numerically integrating the equations to produce 
a height-time relation, which was then used to fit the data. TI fit 
w as implemented with th e hyperbolic equation (Eq. (4]) derived 
by'Kliem & Torok* ('2006'), with and Pi as the fitting parame- 
ters. Poly fit was a polynomial of the following form: 



h{t) = ho-\- vo(t - to) -h ao(t - tof, 



(7) 



where to corresponds to the time of the first data point. The fitting 
parameters were ho, vq, ao and a. The Poly fit not only can pro- 
vide an empirical profile for the later-stage kinematics, but also 
may be used to distinguish the three models because the pro- 
ponents and/or modellers of TI, BO and CA have all proposed 
such polynomial profile with difi'ere nt a for pa rt of their respec- 
tive CME kinematics (~3 for TI; Sc hriiver et al. (2008), ~2 for 
BO (cf. Eq.El), and 2.5 for CA (cf. Eq.^). 

As explained in Sec. 12.11 and 12. 2[ the profiles for CA, 
CAlate.e and TI fits are the profiles for the early stage of the 
eruption. Nevertheless, all four fitting models were applied to 
all the data sets. The reason for fitting the higher coronal kine- 
matics with CA, CAlate.e and TI profiles was that we wanted 
to test the exclusiveness of these profiles to the initial stage. If 
the kinematics of the later stages can also be well-fitted by these 
profiles, the physical assumptions and conditions imposed to de- 
rive them may not be sufficiently strict. Our strategy of using 
these profiles to verify and discern diff'erent models would thus 
be inappropriate and unreliable. The data from diff'erent instru- 
ments were fitted separately because the full kinematics profile 
of a CME process is too complex to be described by one single 
function. In addition, by fitting diff'erent data sets separately, we 
avoided possible errors resulting from instrumental diff'erences. 

The fitting results of CME06 and CME07 are presented in 
Fig- Eland [71 respectively. In the left panels of both figures, we 
plotted only the results of the lower corona to allow better vis- 
ibility of the kinematics in this region. The symbols along with 
the error bars in the figures represent the results from the obser- 
vational data, and the continuous lines are the results from the 
fitting, as denoted in the figures. 
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4. Results 

4.1. 2006 December 17(CME06) 

The initial part of this CME matches the picture proposed by CA, 
that is, the CME rose before the SXR emission began to increase. 
The velocity profile (cf. Fig. [2]) also resembles the model pro- 
fil e that corr esponds to an intermediate reconnection rate (Fig. 6 
in lLin & Fo rbes, 2000). However, a comparison of the kinemat- 
ics and SXR emission profile appears to indicate that there may 
be breakout reconnections before the formation of the current 
sheet. This is explained as follows: If the reconnections in the 
current sheet were indeed reflected in the re- acceleration of the 
CME, the rising in the SXR emission before this may appear 
to be contradictory to the CA model. Although the dissipation 
of the kinetic energy under an ideal MHD process can result 
in thermal emissions, the energy is much lower than the energy 
of X-ray emissions. An explanation we propose to explain this 
is that there may be reconnections between the rising flux rope 
and the overlying magnetic fields before the current sheet be- 
neath the flux rope was formed. While TI theory also proposed 
an ideal MHD process to initiate the instability, most of the kine- 
matic profiles of the CME leading edge from TI simulations did 
not show a decelerati on phase (see, e .g.,lKliem & Torok, 2006; 
iTorok & Kliemll2007l:ISchriiver et al.L l2008). However, we need 
further examinations to determine whether TI may contribute to 
any part of this CME. 

The height-time curve of FIT can be fitted equally well by 
CA, CAlate.e and Poly expressions, but not TI. Hence, only the 
fitting results of CA, CAlate.e and Poly were shown in Fig. [6j 
The discrepancies among the three fittings only became visible 
in the subsequently derived velocity and acceleration profiles. 
Although all three fittings were consistent in the early part of the 
data, only CAlate.e re-produced the fast rising in the accelera- 
tion seen after -14:40 UT. Since the CA-type expressions match 
the observational data, we can use the CA theory to deduce 
that the eruptive flux rope seen in FIT may be free of current 
sheet and that the current sheet may be formed shortly after the 
CME went out of the FIT FOV. The critical height, Aq, obtained 
from our CA fit was approximately 200 Mm, the corresponding 
Alfven speed, Vao, was between 130 and 250 km s"^ and the ini- 
tial speed was roughly 3 km s"^ The observation showed that the 
foot-point separation of this CME was approximately 550 Mm 
(cf. Fig.[T]). Our fitting-determined critical height, being roughly 
half of the foot-point separation, is consistent with the CA theory 
(cf. Sec. 12.1b . Assuming the plasma density p ~ 10"^kg km"^ 
and using Vao ~ 200 km s"^ we can estimate the ambient mag- 
netic field strength at the critical height as follows: 

Bao = VaoV^^ 

^ 7 (G), (8) 

which is a reasonable value at ~ 200 Mm above the limb. 

The exponent a obtained from our Poly fit, which gave a 
good match with the beginning part of the data, was closer to 2 



rather than 3, which diff'ered from the result by Schriiv er et al 
(2008). He nce, the TI expressions produced by Schriive r et al 



(I2008h and iKHem & TorS^ (2006) both failed to fit our data. 
This Poly fit result, however, cannot be used to reject either BO 
or CA because a was mostly scattered between 2 and 2.5. 

The best-fit profile of our best CAlate.e fit had an exponent 
of 5.3. This leads to the following expression for the stage just 
before the formation of the current sheet: 



h : 



(-Vao 

71 



5.3 



(9) 



Aq and Vao given by the fitting were ~ 27 Mm and 100 km s"\ 
which appear to be inconsistent with the values from the CA 
fit. However, since the exponent in the expression is very dif- 
ferent from the theoretically derived value 0.5 (cf. Eq. [3]), Aq 
and Vao from our CAlate.e may no longer represent the critical 
height and Alfven speed. We found that setting the exponent to 
0.5 failed to fit the data. This indicates that while this CME may 
be driven by CA type instability, the assumptions and approxi- 
mations imposed to derive Eq.[5]were inadequate for this CME. 

TI, CA and CAlate.e expressions all failed to fit the LASCO 
data, which is as expected because LASCO only observed the 
later stage of the CME event. This is a promising indication that 
the expressions for our quantitative examination are appropri- 
ate verification tools. Poly fit was the only one that rendered a 
close match to the height-time data. However, some diff'erences 
between the fitting results and the data become visible in the 
velocity and acceleration profiles. Specifically, Poly fit was in- 
capable of re-producing the dip in the C2 velocity. Hence, Fig. [6] 
shows the result of fitting only the later, linear part of C2. In 
addition, the acceleration profile derived from the Poly fit result 
failed to resemble the profile derived from the data. 

The fitting results of FIT and LASCO both revealed that the 
discrepancies become larger in the higher order derivatives. This 
suggests that the derived properties (e.g., velocity and accelera- 
tions) are more discerning than the directly fitted property (e.g., 
height) in verifying a model expression. 



4.2. 2007 December 31 (CME07) 

Although the initial eruption phase of this CME was occulted, 
as we can infer from the derived initial velocity and acceleration 
profile (cf. Fig.O, the peak acceleration in this phase is likely to 
coincide with the maximum of the de rivative of the SXR profile 
(lNeupertlll968l:Plb"mmer et al.Ll2008l) . By a quick visual inspec- 
tion of Fig. O we estimate that the peak occurred approximately 
between 00:45 and 00:50 UT, which is also when a HXR emis- 
sion was detected. 

The kinematics profile of this CME and the SXR and HXR 
emission profiles are in good agreement with a BO scenario: The 
SXR emissions and flaring activities before the emergence of 
the loop and the HXR signal detected on the top of the emerged 
loop could be evidence of breakout reconnections. The second 
HXR signal, coinciding with the beginning of the re- acceleration 
phase, indicated that the reconnections cut ofl' the drag force 
from the current sheet beneath the CME loop and accelerated the 
CME again. However, because the very beginning rising stage of 
the CME may have been obscured by the magnetic arcades, the 
CME may have begun to rise earlier than the increase of the SXR 
emissions, in which case, the timing would be in agreement with 
the CA theory. Since the reconnection is not a necessary con- 
dition to either initiate or accelerate a CME in a TI picture, we 
cannot use the X-ray emissions and flaring activities to verify the 
TI model. However, as we explained in Sec. l4.1[ most TI simu- 
lations of the CME leading-edge kinematics did not show a drop 
of acceleration in the early phase of the eruption. Based on the 
above discussion and reasoning, we deduce that this CME was 
most likely driven by a BO mechanism. 

The kinematics of CME07 were well re-produced by the 
Poly fit, but not the other three models. TI and CA-type expres- 
sions failed to fit even the height data. Hence, only the results of 
Poly fit were plotted in Fig.[7l The unsuccessful fitting using the 
TI and CA-type expressions is another indication that the CME 
was driven by neither of the mechanisms and/or that the early 
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Fig. 6. The fitting results of CME06. The best results of apply- 
ing diff'erent fitting models (i.e., CA, CAlate.e and poly) to the 
height data of EIT and LASCO C2 and C3 are plotted in dif- 
ferent line styles. Diff'erent symbols represent the observational 
results of diff'erent instruments, as denoted in the figure. 



Fig. 7. The fitting results of CME07. The fitting model is a poly- 
nomial of the form Eq.[71 Diff'erent symbols and continuous lines 
respectively represent the observational and fitting results from 
diff'erent instruments, as denoted in the figure. 



stage was not detected, which are both in agreement with the 
results from our qualitative comparison. 

As can be seen in Fig. [71 despite a few mismatching points 
in the velocity and acceleration, overall, the kinematics, both the 
height and the derived velocity and acceleration, of the observa- 
tion and of the fitting were in good agreement. Hence, we can 
deduce that each stage of this CME can be approximated as a 
polynomial of the form h = ho -\- vq ^ t -\- gq ^ even though 
the kinematics of entire process is too complex to be fitted by 
one single function. The fitting results suggest that the order of 
polynomial (i.e., a) for the height-time data from the three in- 
struments were ~ 3 for EUVI, ~ 2 for corl and ~ 2.65 for cor2. 

5. Conclusions 

One of the main goals of this study was to investigate the possi- 
ble driving mechanisms behind coronal mass ejections. We have 
analyzed, both qualitatively and quantitatively, two CME events, 
and compared them with three CME models for this purpose. 

The first CME was observed on 2006 December 17 by EIT 
and LASCO, and the associated X-ray emissions were recorded 
by GOES and RHESSI. This CME lasted for over eight hours, 
and the SXR emission exhibited two rising phases. A compari- 
son of the SXR emission profile and the kinematics of the CME 
showed that the CME was launched earlier than the rising of the 
SXR emissions, reached an acceleration peak during the second 
rising phase of the SXR emission. The CME was first deceler- 
ated after the eruption, and then re- accelerated to another accel- 
eration peak of -60 m s"^, at which time the SXR emissions 



reached the maximum and a flare was detected. The CME prop- 
agated away at a small, constant acceleration at the end. 

Comparing the above information with the three CME mod- 
els discussed in this paper, CA, TI and BO, we propose that this 
CME can be best described by a combination of CA and BO 
models. That is, the CME was initiated by an ideal MHD pro- 
cess (CA), as hinted by the absence of SXR emissions during the 
initial stage. The rising and expanding of the flux rope against 
the ambient magnetic fields then led to reconnections with these 
fields (BO). As it continued to rise, a current sheet was formed 
underneath, which exerted a dragging force to slow down the 
rising motion, and, lastly, the current sheet was cut off' by recon- 
nections and the CME was re- accelerated. Our quantitative ex- 
aminations also revealed that the CA-type expressions gave the 
best fit to the EIT data. From the fitting results, we deduced that 
the critical height at which the CME erupted was approximately 
200 Mm and the Alfven speed at this height was around 130- 
250 km s"\ which leads to a magnetic field of about 7 Gauss. 
They are of reasonable order of magnitude for the low corona. 
Based on the model results by lLin & ForbesI (^000), we can also 
infer from the occurrence of the deceleration of this CME that 
the reconnection rate in the current sheet may correspond to a 
Alfven Mach number between 0.005 and 0.041. 

The second CME was observed on 2007 December 3 1 . This 
CME was formed in a magnetic arcade system. Our analysis in- 
dicated that the initial eruption stage might have been obscured 
by the arcade system. After this un-detected eruption, the accel- 
eration of the CME dropped first and then rose again to a peak 
value of -600 m s"^. The CME remained at this peak acceler- 
ation until its velocity reached a peak value of -1000 km s"^ 
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The CME then decelerated and eventually propagated away at 
a constant speed of -700 km s"^. A comparison between the 
kinematics of the CME and the SXR and HXR emission in- 
formation fr om GOE S, MESSENGER SAX, and RHESSI (cf., 
Fig.S RaftervetaLlJ2010: Krucker et al, 2009a,b) showed that 
the SXR emission began to rise before the emergence of the 
CME and peaked around the same time as the peak of the CME 
velocity, and that the HXR emission due to particle accelerations 
occurred at the time of the CME emergence and at the moment 
when the acceleration of the CME began to rise again. 

These observed features can be well described by a BO 
model as follows: the SXR emission before the CME emerged 
from the magnetic arcades and the high velocity and acceleration 
at the beginning of its emergence can be explained as that the 
breakout reconnections were happening as the CME loop was 
rising through and breaking the overlying arcade. The reconnec- 
tions accelerated the CME, and caused the X-ray emissions. By 
the time the CME broke through the arcades, a trailing current 
sheet may have been formed to act as a drag force, which was 
reflected in the decrease of the acceleration. However, we note 
that some overlying features, such as streamers, could also slow 
down the CME. The subsequent rising of the acceleration (after 
01:00 UT) and the detection of a HXR signal indicated that this 
drag force was being reduced by the reconnections in the current 
sheet. After the current sheet was completely cut ofl', the ambient 
coronal magnetic fields moved in to reform and settle to a lower 
potential state. The last stage of the CME indicated that it was 
being slowed down by the solar wind and/or interstellar medium. 
To be in agreement with a CA model, which proposes an ideal 
MHD process to initiate the eruption, a CME should have a be- 
ginning rising phase during which there was no X-ray emission. 
This is inconsistent with the observations of this CME, which 
detected flares and X-ray emissions from the very beginning of 
the eruption. However, since the initial rising stage of the CME 
was occulted, we cannot completely rule out the possibility that 
an ideal MHD process had taken place. The missing of this ideal 
MHD process in our data was also reflected in our unsuccessful 
fitting using CA-type expressions. 

Both events exhibited a re- acceleration phase, during which 
the maximum acceleration coincided with the peak SXR emis- 
sion. If higher SXR emissions indicate greater reconnection 
rates, as explained in Sec. 12. 11 such coincidence can be explained 
as that the faster reconnections, which lead to faster reduction of 
the tethering force, increased the acceleration of the CME. 

In conclusion, we have investigated two CME events. Our re- 
sults showed that the first event, CME06, can be best described 
by a combination of the CA and BO models while the other 
event, CME07, can be well explained by the BO model alone. 
However, it is also possible that CME07 may have been initiated 
by an ideal MHD process. We conclude that TI model is the least 
likely driving mechanism for either of the events mainly because 
our observationally derived kinematic profiles do not match the 
profiles and expressions from the TI model. Specifically, the de- 
creasing of acceleration in the early st age of our CM E events 
is not seen in TI singulations to-date (iKliem & Toroki 120061: 
iTorok & KHemL 120071: ISchriiver et al.L r2008h . Since our data is 
a two-dimensional projection of the actual motion, we acknowl- 
edge the possibility that this phase be a projection eff'ect. Hence, 
to improve the accuracy of our analysis, we plan to utilize stereo- 
scopic data to obtain an actual, three-dimensional kinematics of 
CMEs in our future work. At last, to improve the current CME 
models, our study suggests that the CA model include the re- 
connections/interactions between the CME and its ambient mag- 
netic fields, BO model consider the possibility of an ideal MHD 



process in the early stage, and that TI model provide possible 
mechanisms for a drag force during the eruptive phase. 
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